Enhanced triplet superconductivity in noncentrosymmetric systems 



Takehito Yokoyama, Seiichiro Onari, and Yukio Tanaka 
Department of Applied Physics, Nagoya University, Nagoya, ^6^-8603, Japan 
and CREST, Japan Science and Technology Corporation (JST) Nagoya, ^6^-8603, Japan 



(Dated: February 6, 2008) 

We study pairing symmetry of noncentrosymmetric superconductors based on the extended Hub- 
bard model on square lattice near half-filling, using the random phase approximation. We show that 
d + /-wave pairing is favored and the triplet /-wave state is enhanced by Rashba type spin-orbit 
coupling originating from the broken inversion symmetry. The enhanced triplet superconductivity 
stems from the increase of the effective interaction for the triplet pairing and the reduction of the 
spin susceptibility caused by the Rashba type spin-orbit coupling which lead to the increase of the 
triplet component and the destruction of the singlet one, respectively. 



Recent discovery of heavy fermion superconductor 
CePtaSi has opened up a new field of the study of 
superconductivity. 1 This is because this material does 
not have inversion center, which has stimulated further 
studies^^i^^^^ Because of the broken 
inversion symmetry, Rashba type spin-orbit coupling 16 
(RSOC) is induced, and hence different parities, spin- 
singlet pairing and spin-triplet pairing, can be mixed in 
superconducting stated From a lot of experimental and 
theoretical studies, it is believed that the most possible 
candidate of superconducting state in CePtaSi is s+p- 
wave pairing £ && 2 && 1£ L I L 1£ In general, d+/-wave pair- 
ing or other mixtures are allowed in noncentrosymmet- 
ric superconductors depending on material parameters. 11 
One important and interesting point about noncen- 
trosymmetric superconductors is how relative magnitude 
of the singlet and triplet components is determined by 
microscopic parameters, which has been unclear up to 
now. 

Spin-orbit coupling causes spin flip scattering and 
hence could influence superconductivity. How the super- 
conductivity is influenced seems to depend on whether 
it is singlet or triplet superconductivity. Recent experi- 
ment of the penetration depth of I^PdaB and I^PtaB 
clearly shows the different behavior of the penetration 
depths, reflecting the difference of the strengths of the 
RSOC. 18 This implies that spin-orbit coupling enhances 
the spin-triplet component of noncentrosymmetric super- 
conductors. However, its physical origin has not yet been 
clarified. 

Motivated by these facts, in this paper, we study pair- 
ing symmetry of noncentrosymmetric superconductors 
based on the extended Hubbard model with the ran- 
dom phase approximation (RPA). As a model calcula- 
tion, we focus on the square lattice near half-filling, which 
is known as a prototype of strongly correlated electron 
systems. We show that d + /-wave pairing is favored and 
that the triplet /-wave state is enhanced by the RSOC. 
Its physical origin is the increase of the effective interac- 
tion for the triplet pairings and the reduction of the spin 



susceptibility which lead to the increase of the triplet 
pairing and the reduction of the singlet one, respectively. 
This could explain the finding in Ref^. 

Let us start from explaining our model. We consider 
square lattice without inversion center. The extended 
Hubbard model including the RSOC can be written as 
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with V(k) = 2V (cos k x + cos k y ). Here, k represents two 
dimensional vector. We set lattice constant to be unity. 
The first term is the dispersion relation and the second 
term represents the RSOC with coupling constant A. The 
third and forth terms represent on-site and nearest neigh- 
bor electron-electron repulsions, respectively. 

Then, the bare Green's functions have the following 
form in the 2x2 spin space: 
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with Matsubara frequency uo n . 

The linearized Eliashberg's equations with RPA in the 
weak coupling approximation are described as 
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with s =|,| and inverse temperature (3. Here, xs and 
Xc are spin and charge susceptibilities at uo n = 0, re- 
spectively, which are obtained by xs = X^ ~ Rex^ 
and xc — X^ + Rex^- Note that x^ = X^ an d 
X^ = (x^) are satisfied, x^ an d X^ are given by 

/x TT (fc)\_ if xi + v(k)(xl-\xi\ 2 ) \ 
U n M)~A{xi-(u + v(k))(xi-\xi\ 2 )) ,{) 

A = 1 + 2^ Xl -U(U + V (A;)) (x? - |X2| 2 ) + 2(17 + 
^(fc))Rex2, 

Xi (k) = -jjy ^2 G + (k + q,ioJ n )G+ (q,iu n ), (10) 
X2(k) = --^ ^ G n (fc + g, zo; n ) G|t ^n). (11) 
xL T d and xjad are defined as 

B = (l- Uxif ~ \Ux2t Notice that X \ld = x\L ™d 
xJid = (xL) are satisfied. 

In the RPA, we take into account the contributions 
from bubble and ladder types of diagrams.— In the above, 
X^ and x^ stem from the bubble type of diagrams, while 
xjad an d Xiad or igi na ^ e from the ladder type of diagrams. 
Here, we ignore the nearest neighbor electron-electron 
repulsion in the ladder diagram for simplicity. However, 
we can grasp the essence of the physics by the RPA . 2Q i 21 

By solving the Eliashberg's equations, we can obtain 
the gap functions. We define singlet component of the 
pair potential and triplet one with S z = for a later 
convenience as 

A s = (A n -A iT )/2, (13) 
A t = (A u + A iT )/2. (14) 
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FIG. 1: (color) Real and imaginary parts of gap function A^. 
We take A = 0.01 and A = 0.5 in the upper and lower figures, 
respectively. Solid lines represent Fermi surfaces. Arrows 
indicate typical scattering processes. 



Let us first discuss the pairing symmetries. We put 
t = 1 and set parameters as T = 0.01, N = 256 x 256, U = 
1.6, n = 0.9, and V = 0.3, while we change the strength 
of RSOC A. Here, T, A/", and n denote the temperature, k- 
point meshes and the band filling, respectively. For these 
parameters, the eigenvalue A is around A = 0.3 ~ 0.7. 
Gap functions are normalized by the maximum value of 
|A|J. Note that An = — A*^ and A s G R are satisfied 
due to the time reversal symmetry. We have found that 
the magnitude of the triplet component with S z = 0, 
A t , is negligibly small. Then the d- vector has no z- 
component. This is consistent with the prediction that 
the d- vector tends to be parallel to (sin k y , — sin k x , 0). 4 
We have also confirmed that the obtained results do not 
change qualitatively for other set of parameters. 

Figure [1] depicts real and imaginary parts of gap func- 
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FIG. 2: (color) Real part of singlet gap function A s . We 
take A = 0.01 and A = 0.5 in the upper and lower figures, 
respectively. Solid lines represent Fermi surfaces. 



tion Ajf. As shown in this figure, it has a /-wave sym- 
metry. We take A = 0.01 and A = 0.5 in the upper and 
lower figures, respectively. We can see that the magni- 
tude of the gap function is enhanced with the increase of 
RSOC. 

We show the real part of singlet gap function A s in 
Fig. El We also take A = 0.01 and A = 0.5 in the upper 
and lower figures, respectively. As can be seen from this 
figure, it has a <i-wave symmetry, which is consistent with 
the previous works j 2Q i 21 We can find that the symmetry 
of the gap function is independent of the RSOC. 

The appearance of /-wave symmetry in at A = 
0.01 can be understood by the structures of the spin and 
charge susceptibilities (xs and xc)- They have peaks 
near (±7r, ±7r) as shown in Fig. [3l According to Eq.([7j), 
the gap functions tend to have the same sign during 
the scattering process. Therefore, /-wave symmetry is 
favored. We show typical scattering processes toward 
(7r, 7r) by yellow arrows in Fig. [TJ A rather complicated 
structure of the real part of A^ at A = 0.5 in Fig. [1] 
stems from the fact that the off-diagonal components of 
the Green's functions disappear at the van Hove singu- 
larities (see Eq.Q). Since the contributions from the di- 
agonal (off-diagonal) components are dominant near (far 
away from) the van Hove singularities, we can expect line 
nodes in the intermediate regions. 

To understand the origin of the enhancement of the 
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FIG. 3: (color online) Spin and charge susceptibilities (xs 
and Xc)- We set A = 0.01 and A = 0.5 in the upper and lower 
figures, respectively. 



triplet pairing, we study xs and xc-> and the real part of 
Xiad' wn i cn are plotted in Fig. [3] and Fig. [H respectively. 
Clearly, xs is reduced by RSOC while xc is almost in- 
dependent of it. This can be intuitively interpreted as 
follows. Spin-orbit coupling causes spin flip process and 
hence breaks magnetic fluctuation. On the other hand, 
spin flip scattering does not affect charge fluctuation. 
Therefore, xs depends on the RSOC while xc is almost 
independent of it. Since the positions of the peaks in 
Xs and xc are almost the same, they compete with each 
other (see Eq.(|5J)). It is known that in such a situation, 
the decrease of xs leads to the reduction of singlet pair- 
ings and hence triplet pairings could dominated Note 
that xlad i s a ^ so recm ced by the RSOC. Additionally as 
shown in Fig. 0J the real part of xjad^ which contributes 
to the effective interaction for the triplet pairings, is en- 
hanced by the RSOC. This results in the enhancement 
of the triplet component. Note that the imaginary part 
of xJad i s negligibly small. These are the reasons of the 
enhancement of the triplet pairing by the RSOC. 

Finally, let us discuss the relative magnitude of 
max | A^ | and max|A s | in detail. In Fig. [5j we plot 
max |A|f | / max |A S | as a function of the strength of 
RSOC A. It increases monotonically with A. For a suffi- 
ciently large magnitude of A, the magnitude of the triplet 
pairing could be comparable with the single one. 

In conclusion, we studied pairing symmetry of non- 
centrosymmetric superconductors, where we used the ex- 
tended Hubbard model on square lattice near half-filling 
with the RPA. We found that d + /-wave pairing is fa- 
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FIG. 4: (color online) Real part of xlad- We se ^ ^ = an< ^ 
A = 0.5 in the upper and lower figures, respectively. 




FIG. 5: Relative magnitude of max |Ajj| and max |A S | as a 
function of the strength of RSOC A. 



vored and that the triplet /-wave state is enhanced by 
the RSOC. This stems from the increase of the effective 
interaction for the triplet pairing and the reduction of the 
spin susceptibility by the RSOC, resulting in the increase 
of the triplet pairing and the destruction of the singlet 
one, respectively. 

Our results are consistent with the recent ex- 
periment of the penetration depth of I^PdaB and 
Li 2 Pt3B which indicates that spin-orbit coupling en- 
hances the spin-triplet component of noncentrosymmet- 
ric superconductors^ 
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